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The scattering theory of quantum transport relates transport properties of disordered mesoscopic 
conductors to their transfer matrix T . We introduce a novel approach to the statistics of transport 
quantities which expresses the probability distribution of T as a path integral. The path integal 
is derived for a model of conductors with broken time reversal invariance in arbitrary dimensions. 
It is applied to the Dorokhov-Mello-Pereyra-Kumar (DMPK) equation which describes quasi-one- 
dimensional wires. We use the equivalent channel model whose probability distribution for the 
eigenvalues of TT^ is equivalent to the DMPK equation independent of the values of the forward 
scattering mean free paths. We find that infinitely strong forward scattering corresponds to diffusion 
on the coset space of the transfer matrix group. It is shown that the saddle point of the path integral 
corresponds to ballistic conductors with large conductances. We solve the saddle point equation and 
recover random matrix theory from the saddle point approximation to the path integral. 

PACS numbers: 72.10.Bg, 05.60.+W, 72.15.Rn, 73.50.Bk 



I. INTRODUCTION 

Advances in micro-fabrication technology led to the re- 
alization of mesoscopic electronic devices. In such de- 
vices the mean free path for inelastic electron scattering 
exceeds the dimension of the device. As a consequence 
phase coherence is maintained which leads to quantum 
interference effects like universal conductance fluctua- 
tions, persistent currents and Aharanov-Bohm oscilla- 
tions in rings, or weak localization!!. The phase coher- 
ence has also serious theoretical implications. It causes 
large conductance fluctuations which are related to the 
problem of high gradient operators in the field ptheoretic 
description of the metal insulator transitiona cl These 
fluctuations manifest themselves already in the metallic 
regime as logarithmic normal tails of the conductance 
probability distribution. As the critical regime is ap- 
proached the conductance probability distribution be- 
comes increasingly broader until it reaches a logarithmic 
normal form in the insulating regimeO. 

A common approach to transport quantities of meso- 
scopic conductors is the scattering theory of quantum 
transportLHa. It models the conductor by a disordered re- 
gion which is connected to a number of ideal leads which 
support propagating wave modes. The number of leads 
corresponds to the number of measurement terminals. 
Here only two terminal geometries will be considered. 
The scattering matrix relates the amplitudes Ik , I' k of the 
incoming with the amplitudes Ok,0' k (k = 1, . . . , N) of 
the scattered propagating wave modes at the Fermi en- 
ergy, 



S = 



O 

a 



(i) 



r f 
t r' 



(2) 



t and r are the transmission and reflection matrices for 
incident waves from the left, and t 1 and r' are the trans- 
mission and reflection matrices for incident waves from 
the right. The dimensionless two-probe conductance 
g = G/(e 2 /h) in terms of the transmission eigenvalues 
T fc of tit i s 



(3) 



k=l 



where 



There are three universality classes which correspond to 
different physical situations. Conductors with time re- 
versal invariance lie in the orthogonal universality class. 
The unitary universality class corresponds to conductors 
in which the time reversal symmetry is broken, e.g. by a 
magnetic field. Conductors with spin-flip scattering pro- 
cesses but no time reversal symmetry breaking fall into 
the symplectic universality class. 

Recently the quasi-one-dimensional wire has attracted 
considerable attention. The width of a quasi-one- 
dimensional wire is of the order of the mean free path 
for elastic electron scattering so that transverse diffusion 
can be neglected and the cross section of the wire be- 
comes structureless. Interesting non-perturbative results 
which are valid-dpe all wire lengths have been obtained 
for this systerrollj. Furthermore it has been the ideal 
playground for new ideas in the field of quantum trans- 
port. 

One of these ideas is the Fokker-Planck (FP) ap- 
proach to quasi-one-dimensional wires. The FP equation 
which describes the probability distribution for the trans- 
mission eigenvalues is known as the Dorokhov-Mello- 
Pereyra-Kumar (DMPK) equation. It has been derived 
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by a number of author sli3~tLZI who started from various 
different models. Its form is 



dp(s;{T k }) 



ds 



dp 



(3p 



dn({T k }) 



(4) 



where 



^({r fc }) =-^ln|(coshr fc -coshr z )/2| 

k<l 

-l//3^1n|sinhr fe |, 



(5) 



7 = /3N + 2-/3, and coshT fc = (2-T k )/T k . The values of 
(3 are 1,2, and 4 for the orthogonal, unitary and symplcc- 
tic universality class respectively. The DMPK equation. 

ly in the past few yearsE3~L3. 



has been studied inte, 
Beenakker and Rejaeir 8 ' 



discovered that the variation 



p(s; {T k }) = ex P |-|n({r fc })| ^(a; {r fc }) (6) 

of the Sutherland transformationS which is known to 
solve theJirownian motion model for the circular unitary 
ensembleEil, works as well for the DMPK equation. Af- 
ter this transformation ip(s; {F k }) obeys a Schrodinger 
equation for TV non-interacting particles. As a conse- 
quence the exact form of p(s; {T^}) could be determined. 
This solution has been the basis for Frahm's exact calcu- 
lation of the one- and two-point correlation functions of 
the transmission eigenvalues^. 

In this paper we present a novel approach to the scat- 
tering theory of quantum transport which expresses the 
probability distribution of the transfer matrix as a path 
integral. Our motivation has been the belief that the 
path integral technique can be developped into a tool 
which is more powerful than the FP approach when it 
comes to the description of higher-dimensional conduc- 
tors. 



II. SCATTERING MODEL 

We use the transfer matrix T instead of the S-matrix 
to model the scattering properties of the disordered con- 
ductor. The transfer matrix relates the scattering ampli- 
tudes in the left lead with the scattering amplitudes in 
the right lead 



O' 
V 



(7) 



It has the advantage that it obeys the multiplication law 
T(L + SL,0) = T(L + 5L,L)T{L,0) (8) 
which leads to the simple Langevin equation 



T(x)= \ ' = e(x)T(x,0) 



dx 



z iL (x) e i2 (x) 
z 21 {x) e 22 (x) 



T(x,0) (9) 



for the stochastic evolution of the transfer matrix. The 
disorder is generated by the multiplicative noise e. 

In this paper we consider only conductors in the uni- 
tary universality class. Then, T obeys the symmetry 
constraint 

£ z Tt£ z T = l (10) 
which ensures flux conservation, where 



S 2 = 



1 
1 



(11) 



A convenient parametriaation of the transfer matrix is 
the polar decompositioiJl§tB 

T 



m 

u 3 



cosh(r/2) sinh(r/2) 
sinh(r/2) cosh(r/2) 



U2 
«4 



(12) 

where T is a real, diagonal N x N matrix and Ui (i — 
1, 2, 3 and 4) are unitary N x N matrices. 

The relation ( |To| ) implies that S 2 e^S z +e = leading 
to the symmetries 



jilt = 
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e 1 = e 

for the noise. The stochastics property 



(13) 
e could be 

derived from a microsmpic HamiltonianOE.]. Here, we 
adopt a simple mode£jl£j which assumes Gaussian white 
noise such that 



(s k }(x)s k } l r(x')) 

(elUx)e k %f(x')) 
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S kk '5 w 6(x - x'), 



(14) 



and all other independent second moments are zero. The 
mean free paths Zjy, l'^, and l\p l' h kl for forward and back- 
ward scattering, respectively, are defined by the limits of 
the disorder averages 
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for a short piece of conductor with length SL. Note that 
the symmetries (13) imply the relation l h kl = l'f k . 

We want a path integral representation of the stochas- 
tic process (^) in terms of the transfer matrix T. The 
derivation technique which is most suited for that pur- 
pose derives the path integral directly from the Langevin 
equation (see chapter 4 in Ref. |3q ). The symmetry 
constraints ([To]) on T will be taken into account by 6- 
functions which leads naturally to the invariant measure 
of the transfer matrix group as the path integration mea- 
sure. We illustrate the essential ideas of the derivation 
technique with the simple example of diffusion on a circle 
before we deal with the transfer matrix. 



III. DIFFUSION ON THE CIRCLE AS A SIMPLE 
EXAMPLE 

Let the angle ip determine the position on a circle. The 
analogue of the Langevin equation (||) is 



u = ™ - e{t)u{t) 



(16) 



where u = exp(iip). The symmetry e* = — e implies 
d(uu*)/dt = which ensures that u remains a phase. 
Choosing Gaussian white noise for the imaginary part of 
e such that 



(e(t)) = 0, 
(e(t)e(t')*) =2D5(t 



n 



leads to the FP equation 



dt 



D- 



{dp? 



(17) 



(18) 



which describes diffusion on the circle. 

The probability distribution of u can be formally ex- 
pressed as 



p(t;u) = (5(u-u(t))) 



(19) 



where u = + iu^ 2 \ S(u) = S(u^)S(u^ 2S> ), and u(t) is 
the value of u which is acquired at time t for a certain 
realization of the noise and the initial value u(0) = uq. 
The brackets (• • •) denote the average over all possible 
noise configurations. The path integral representation is 
derived by inserting a product of 5-functions 



p(t;u) = ( / ]J du(t')S{u(t') - u(t'))S(u{t) - 



t'=0 



(20) 



where du = du^du^. The (5-function S(u(t') - u(t')) 
restricts the value of u(t') to u{t'). Since u(i') is not 
explicitly known we enforce this constraint implicitly by 



the relation ii(t)u 1 (i) — e(t) — which follows from the 
Langevin equation (|16[). That leads to 



p(t;u) = ( / II du{t')\detA\5{u{t')u-\t') - e{t')) 

x6(u(t)-u)\ (21) 



t'=0 



where the operator A is defined by the functional deriva- 
tive 

^ _ S(u(t)u-\t) - s(t))^ 
Ayl ' >- 6uU')(t>) ■ 1 ' 

The average over the Gaussian probability measure 
P[e] n de{x) 4exp(-f dx £{x) ** {x) ) f[ de(x), 

x=0 [ Jo ) x=0 

(23) 
(24) 



where 



de = de w de (2) S(e + e*), 



yields 



j(t;«) = ^ f \[du{t')5{u{t')u- 1 {t')+u*{t') 

J +>—r\ 



xu-' *{t'))\ detA\ exp{-S}, (25) 



where 



S = 4D J ^(Ott- 1 ^)^" 1 ^))' (26) 

and the path summation includes all paths which start 
at Mo an d end at u. 

The property that w^w^ 1 (t) + w* (i)w _1 * (t) = 
+ u*(t)u _1 *(£) if w(t) = u(t)v(t) and v(t) 
is a phase, suggests that Ilt'=o du(t')5(ii(t')u~ l (t') + 
■U*(t')u _1 *(t')) is proportional to n*'=o dfi(u(t')) where 
d/j,(u) is the invariant measure on U(l). This becomes 
explicit if the (5-function is introduced via an auxiliary 
field rc(t') 

p(t;u) = J Y[du(t')dK(t')\dctA\cxp{-S}, (27) 



t'=o 



where 



S = S + i / dt' K{t'){u{t')u- l (t') + u*{t')u~ l *{t')) 
Jo 



r* d 

= S + i dt'K(t')-^ln(u(t')u*(t')). 
jo dt 

Partial integration yields 



(28) 
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S = S + i dt'\(t')]n{u{t')u*(t')), (29) 
Jo 

where X(t) = —k(t). The Jacobian of the transforma- 
tion A(t) = —k(t) is an irrelevant constant which can be 
incorporated into the normalization factor. Hence 

P (t; u) = TV -1 f Y[ dn(u(t'))\ det A\ exp{-S} (30) 

J t'=Q 

since duS(hi(uu*)) = duS(uu* — ll which is proportional 
to the invariant measure d/i(u)l2j. The restriction to 
uu* = 1 in the invariant measure simplifies the action 



(31) 



To calculate dot A we evaluate Eq. (|22|) which gives 

An(t,t') = (a(t,t')+a*(t,t'))/2, 

A 12 (t,t') = i(a(t,t')-a*(t,t'))/2, 

A 21 (t,t') = -i(a(t,t')-a*(t,lf))/2, 

A 22 (t, t') = (a(t, t') + a*(t, 0)/2, (32) 

where 

a{t,t')=u- 1 {t) (^j t 5(t-t')-5{t-t')u{t)u- 1 {t)^ . 

(33) 

The decomposition A = BCT> into a product of three 
operators 



[B] (t,0 

[C] (M') 

[£](*, 



i 

71 



<5(t-t') -i8(t-t') 
-iS(t-t') 5{t-t') 



a(t,t') 
a*(M') 



J_fS(t- t') iS(t - t') 
y/2\m- If) 5(t-t>) 



(34) 



implies that det A = det C = det S det a* since det B — 
det D = 1. The operator a can be as well factorized into 

a = 0,10,20,3 where 

ai (*,*') = u -1 (»5(i-*') 



a 2 (t,t') = -S(t-t') 

a 3 (t, t') = S(t - t') - 0(t - t'W^it') 



(35) 



The determinant of SiS^ is one since the 5-function in the 
path integration measure enforces that u(t)u*(t) = 1. 
The determinant of d 2 is an irrelevant constant which 
contributes only to the normalization. Using det = 
exptrln and ln(l + x) — l) k+1 x k /k to evaluate 

det S3 yields 



det S3 = exp 



(36) 



The higher order terms which are indicated by the dots 
vanish due to products of ^-functions. The quantity 9(0) 
is not defined yet which can be traced back to the mul- 
tiplicative noise in the Langevin equation (fHf). The cor- 
rect choice is 9(0) = 1/2 (see the discussion in chap- 
ter 4 of Ref. Eft). Here this choice does not matter 
since u(t')u- 1 (t r ) + ii*(t')u- 1 *(t') = which implies that 
det a 3 det S3 — 1 , leading to the final form 

p(t; u) = Af- 1 f Y[ dn(u(t')) exp{-S*} (37) 

of the path integral representation of the stochastic pro- 
cess (jig). 



IV. THE PATH INTEGRAL FOR THE 
TRANSFER MATRIX 

The analogue of Eq. (|2(]) for the transfer matrix is 
p(L-T) = [( [ f[ dT(x)S(T(x) - f(x)) 

x5(T(L)-T)\, (38) 



x=0 



where 



dT^l[ dT M dT ^ 

k,l 

S(T-T)^]Jd(T^-f t 



kj 



xHT^-f^). 



(39) 



Enforcing T(x) by T(x)T- 1 (x) - e(x) = which follows 
from the Langevin equation (|{j) yields 

p(L;T)= jlj Y[dT(x)\detA\S(T(x)T- 1 (x) - e(x)) 

xS(T(L)-T)\, (40) 



where the operator A is defined by the functional deriva- 
tive 



, s _ S[T(x)T- 1 (x)-e(x)] U) 



kl 



(41) 



Performing the average over the Gaussian probability 
measure 
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P[e] J] de(x) = 1 ^p[-\J q L dx(l f kl e k }(x)e k }*(x) 
+l'l l ef l {x)ef l *{x)+l b H Ef l {x) £ f l *{x) 
+l%el}(x)ei}*(x))\f[de(x), 

' x=0 

(42) 

where 

de^ J] d4 (1) 4 (2) *s(*), 

s s (e) = nW( e « + 4 1 *) (1) )<(4? +4 1 *) (2) ) 

x<5((e^+ £ ^*) (1) )<(4 2 +4 2 *) (2) )} 

n<(4 i fc ) (i) )<ra (i) ) 
fe 

(43) 

yields 

p(L;T)=Af- 1 / JJdT^^^T- 1 ^) 

«* n 



x=0 



x|det„4|exp{-S}, 

(44) 



where 



-ll22 
kl 



+r^[tr- 1 ] 2 fe [[rr- 1 ]^*}. 

(45) 

By analogy with the previous section we expect 
that Y\^ dT(x)S s (f(x)T- 1 (x)) is proportional to 

Yl^o dfj,(T (x)) where d/i(T) is the invariant measure of 
the transfer matrix group. This will be proven in ap- 
pendix [A]. The form of the invariant measure in terms of 
the polar coordinates (O) is 



k<l 



dfi(T) = JJ(coshr fe - coshr ; ) 2 ]^[sinhr fc dr fc 

4 

Y[ dfi(u k ), (46) 



k=l 



whero-6^t(ttfc) is the the invariant measure on the unitary 
grouped. 



We proceed with the calculation of det„4. Using 
d/dT^ = d/dTu + d/dTfr, d/dT$ = i{d/dT kl - 
dJ&Tfa), and dT^ x ldT k , v = -T^ltf to evaluate Eq. 
© yields 

[A] 11 = (A + A*)/2, 

[A] 12 = i{A- A*)/2, 

[A] 21 = -i(A- A*)/2, 

[A] 22 = (A + A)/2, (47) 

where 

Akl,k'l'(x,x') = hrnT M 1 {x)(^-5(x - x')S m k'S n l' ~ 

S{x - a/Ht^T-^U^A (48) 



The decomposition A = BCT> into a product of three 
operators 



B = 



C = 











( A 




V o 


*)• 




' i *i 


71 ( 


v <i i 



(49) 



where [l]w,fc'/'(af,af') = <5(x - x')S kk >5w, implies that 
det A = det (7 = det A det A* since detS = det£> = 1. 
The operator A can be as well factorized into A = 
A1A2A3 where 

Ai.H.k'V&x') = [l(g> {T-Yi^kk'M'Six - x') 

A 2 -ki.k'i'{x,x') = -^-S(x - x')6 k k'Sw 
ax 

A 3 -ki,k'i' {x, x') = [6(x-x')l®l- 9(x - x 1 ) 

xr^Or-^^O^ilfei.fe'i'- (50) 

The product det A\ det is one since the determi- 
nant of the transfer matrix is a phase. The determi- 
nant of Ai is an irrelevant constant which contributes 
only to the normalization. Using det = exptrln and 
ln(l +x) = Efcli(- 1 ) fe+la; ' £ / fc to evaluate det A 3 yields 

det A 3 = exp|-7V6»(0)^ dx t^T^T" 1 ^)) j (51) 

The symmetries of the transfer matrix imply that 
tr(T(a;)T- 1 (x) + (T^T' 1 ^))*) = which gives 
det A 3 det A 3 = 1 . 

That leads to the final form 

p(L; T) = AT 1 / H d^{T(x)) cxp{-S} (52) 

J x=0 

of the path integral, where S is the action of Eq. (E5I) . 
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V. THE DMPK EQUATION 

We formulate the DMPK equation in terms of diffu- 
sion on the coset space of the transfer matrix group as 
has been done by HuffmannEl In our context that can 
be achieved with the equivalent channel model (ECM). 
This model has been introduced by Mfilk) and Tomsovic 
for the orthogonal universality classnZrEJ. They showed 
that it is equivalent to the DMPK equation with (3=1, 
in the sense that the joint probability distributions for T 
of both models are identical. The ECM for the unitary 
universality class is just the model ( |l4| ) with backscat- 
tering mean free paths of the form 



l b 

III II 



1 

IN 



(53) 



and arbitary forward scattering mean free paths. It is 
equivalent to the DMPK equation with (3 = 2 in the 
same sense. The difference between the DMPK equation 
and the ECM is that the unitary matrices need not be 
isotropically distributed and that there can be correla- 
tions between them and T. 

We choose forward scattering to be infinitely strong so 
that the mean free paths l( nn and are zero. Then, 
the action ( |45| ) simplifies 

S=™£dx ^{[TT-^^TT- 1 ]^ 

+ [TT- 1 ] 21 ([TT- 1 ] 21 ) 1 }. 

(54) 

Using that TT 1 = -TT Ll and the symmetries of 
TT -1 one can simplify further 

S =^l J L dx tr{ (TT- 1 + (TT- 1 ) ! ) 2 } 
J dx trj 2 TT- 1 (TT- 1 )' i 



Nl 



Nl 



dx 



tr (MM- 1 ), 



(55) 



where M = T^T which does not depend on u x and u 3 
anymore. The infinite strong forward scattering immedi- 
ately randomizes the probability distribution of U\ and 
M3 so that they become isotropically distributed. Note 
that the space which is formed by the matrices M is iso- 
morphic to the coset space of the transfer matrix group. 
The path integral describes diffusion on the coset space 
since the action is the classical action for free motion on 
this spacer 

Introducing the dimensionless length s = x/(Nl) yields 



S 



1 



l/9c 



ds tr( MM- 1 



(56) 



where the dot now stands for the derivative with re- 
spect to s rand, g c i = Nl/L is the classical (bare) 
conductanceE-St3 in units of e 2 /h. Hence, large con- 
ductances correspond to the 'short time' regime of the 
path integral which justifies a saddle point approach 
for good conductors. The variation M(s) + SM(s) = 
ST^(s)M(s)ST(s), where ST = 1 + e and e obeys the 
symmetries (|l^) leads to the saddlepoint equation 



= 5S' 



l/9ci 



ds tr^ M + Me) M- 1 

-MieM- 1 + M- 1 e t )). (57) 



One can verify easily that M sp (s) = exp{sX} is the so- 
lution for a path which starts at M(0) = 1 and ends at 
M = cxp{X/g c i}. Evaluation of the saddle point action 
yields the transfer matrix probability measure in saddle 
point approximation 



p(i;T)dMT)«n ex p{-^ r '}^( T ) 

= ^(coshl^ - coshT/) 2 JJexp 

k<l k 
4 

J^sinhrfcdrfc dfi(u k ). 



Nl r2 



(58) 



k=l 



This is just the random matrix theory probability distri- 
butionrBieawire which has been proposed for the transfer 
matrixc3c2l'Ei Since it is known that random transfer- 
matrix theory describes the stochastic properties of bal- 
listic conductorsEJ we conclude that the saddle point of 
the path integral correctly describes the ballistic regime 
of the conductor. 



VI. CONCLUSION 



In summary we have presented a path integral ap- 
proach to the stochastic properties of mesoscopic dis- 
ordered conductors. Its application to quasi-one- 
dimensional wires in the ballistic regime led to the ran- 
dom transfer matrix theory probability distribution. We 
believe that known results for the quasi-one-dimcnsional 
wire could be recovered by a systematic perturbation ex- 
pansion in powers of l/g c i- At the moment it is not clear 
to us wether the 'short time regime' of the path integral 
in higher dimensions corresponds as well to conductors 
with large conductances. That still has to be clarified. 
The further development of the path integral technique 
also remains to be done. 
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APPENDIX A: THE INVARIANT MEASURE OF 
THE TRANSFER MATRIX GROUP 

The invariant measure on the transfer matrix group 
does not change under multiplication with a fixed trans- 
fer matrix Tq from the left or the right 



d/*(T) = d(i(T T) = dfi(TT ). 



(Al) 



In this appendix we prove the claim of sect. IV 
that T\t=o dT ( x ) S s(T{x)T- 1 (x)) is proportional to 

U L x =oMT(x)). 

Since the inverse of T in 6s cannot be handled as eas- 
ily as it -1 in the example of diffusion on the circle, we 
show first that 63(e) oc 6s('E z T^'E z £T) up to a Jaco- 
bian. This will allow to replace TT _1 in the argument 
of 5 5 by£ z Tt£ z t. 

Writing the (^-function in terms of its Fourier represen- 
tation yields 



Ss(e) = 



1 



(211) 



4jV 2 



G?/«exp<j -tr[/c(e + £ z e^£ z )] 



where 



K = 



(A2) 



(A3) 



K 11 t = K 11 
« 22 t =K 22 
K 12 t = -K 2 \ 



(A4) 



and 



dK = JJ dnl] m dKu (2) d/s« (2) 



k<l 



.11(1) ,22(1) 



\\dn 

k 

Then the linear transformation 



kk dK k l^'Y[dK 

k k,l 



12(1) ,12 (2) 
kl aK kl 



(A5) 



e' = £ z T f £ z eT 
of e can be absorbed into k, 

Ss(e') = ' 



(A6) 



(2n)4 jv2 J dKexp|-tr[K'(e + S,et Sj( )] 

(A7) 

where 

k' = TkH z T^T, z . (A8) 

Since k! has the same symmetries as K it follows that 

6 s (b')=6 s (b)/\J(T)\, (A9) 

where J(T) is the Jacobian of the linear transformation 
(|AS|). Hence, replacement of the argument TT~ Y in 6s 
by S 2 T^S 2 T via the linear transformation ([A 6]) yields 



\\6 s {t(x)T-\x)) « 

L 

\\dK(x)\J(T(x))\ 



x=0 



x=0 



x exp^ 1 I dxtr 

lx=0 



dx ' 



(A10) 



Partial integration and using that £ z Tt£ z T = 1 at the 
endpoints gives 



x=0 



Y[6 s (T(x)T-i(x)^ 

L 

l[dn(x)\J(T(x))\ 



expj-i J dxtv[k(Y, z T^Y, z T - 

(All) 



The Jacobian of the transformation k = — k is a con- 
stant. Hence 



x=0 



l[6 s (T(x)T- 1 (x)) ex 

L 

J] |^(T(x))|<5 s (S z Tt( a; )S z T(x) - 1) 



x=0 



(A12) 



In order to calculate J(T) we introduce the (4A^ 2 )- vec- 
tor notation 

K T = (Ku, ■ ■ ■ , K12N, «21, ■ ■ ■ , K 2 N2n) (A13) 



7 



of the matrix k. Then %' = (T ® (£ 2 Tt£ 2 ) T ) k. There 
is a complex matrix E such that k = EiZi n d, where K- m d 
contains the AN 2 real and imaginary parts of the inde- 
pendent matrix elements of k. Therefore 



E- 1 (T ® (£ 2 Tt£ z ) T ) En ind . (A14) 



J(T) is the determinant of this linear tr ansfo rmation, 
which is one since the (^-functions in Eq. ( Al2| ) enforces 
S.r+S, to be the inverse of T. That leads to 



L L 



x=0 



JJ*s(E«rt(ar)S,r(a;)-l). 

(A15) 



It remains to be shown that 

d(i(T) = dT5 s (X z Ti-E z T-l) 



(A16) 



has the properties ( jAlp and therefore is the invariant 
measure. 

For multiplication with a transfer matrix T from the 
left the argument of the (5-function does not change which 
leads to 

d(i(T T) = dT\I(T )\6 s (X z T^ z T- l), (A17) 

where I(Tq) is the Jacobian of the linear transformation 
T" = TqT. Expressing this transformation in terms of 
real vectors yields 



f,(i) 

fl(2) 




,(2) 





(A18) 



The Jacobian I(To) is the determinant of the transfor- 
mation matrix which can be decomposed into the product 



1 



-n 



V2 V- Jl 



T <e> i o 

TX <g 



1 



1 il 



V2V a 1 



(A19) 



of three matrices. Since H z T' q Ti z Tq = 1 implies that 
detT detT5 = 1 one finds that I(T ) = 1 and there- 
fore dn(T T) = dfi{T). 

Analogously it can be shown that the Jacobian for the 
multiplication with To from the right is one as well which 
gives 

dfx(TT ) = dT5 s (£ 2 Tt£ 2 (£ 2 Tt£ 2 T- l)T ). (A20) 
As shown above 5s(S 2 TqS 2 £To) = <5g(e). Hence 

d(i{TT ) = dT6 s (X z T^ z T - 1) 

= d^(T) (A21) 
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